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Abstract: We present semi-analytical solutions of the supersymmetric non-universal masses mod- 
els for low tan/? regime. In addition to this, scale and tan/3 dependencies of the soft (mass)^ terms 
are given in the form of numerical solutions. By using the constrained form of the semi-analtic 
results, particular attention is paid on the non-universality assumption of the Higgs mass values 
and their potential measurable effects on the mass spectra of the minimal supersymmetric standard 
model. It is observed that, certain measurables are almost insensitive to the initial mass choices of 
the Higgs fields, like the mass of the light CT^-even Higgs boson. On the other hand, large deviations 
exist on the mass of the remaining physical Higgs bosons signal that the allowed parameter space of 
the model can be probed successfully. For this aim, in addition to the other physical Higgs bosons, 
imprints originating from the heavier chargino {x^), heavy neutralinos (X3, X4) and the light scalar 
tau (fi) are necessary and found to be promising. 
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1. Introduction 

There are a number of motivations for phenomenological studies of the Supersymmetric (SUSY) 
theories among which unification of the gauge couphngs and natural suppression of the radiative 
corrections on the masses of Higgs bosons can be mentioned (see i.e. 0], for a comprehensive list of 
motivations). Among those theories, due to least number of particles, the Minimal Supersymmetric 
Standard Model (MSSM) occupies a special place. In the near future, forthcoming experiments may 
reveal that the incorporation of the Standard Model (SM) into a more effective theory turns out as 
the MSSM. Indeed, if low energy supersymmetry is realized in Nature, phenomenological studies 
related with the MSSM and its variants will be important to unreveal the hidden model. Since it 
has certain problems like the famous /i problem 1^ , flavor problem |^ , and the unknown mechanism 
of the supersymmetric symmetry breaking, studies related with the extensions of the MSSM may 
be expected to shed light on future measurement, especially if nontrivial data inconsistent with the 
minimal model occurs. 

In this work, wc study particle spectrum in the MSSM with non-universal Higgs mass terms 
(NUHM) 1^. We provide most general semi-analytic solutions of evolving terms, in terms of high 
scale boundary conditions, for a low tan (3 value. Additionally, different scale and tan f3 dependencies 
of the soft (mass)^ terms will be presented numerically. Actually, the exploration of solutions to 
the renormalization group equations (RGEs) of a supersymmetric model with NUHM is a subject 
that has been investigated (see e.g. Q and [^), but, the novel feature of our analysis is that 
semi-analtic solutions may facilitate the exploration of the phenomenology of the model (see for 
phenomenology of NUHM). As is well known, weak scale observables and Grand Unified Theory 
(GUT) scale boundaries are connected via RGEs in a complicated manner |^ and they can be 
solved with the help of certain softwares. Taste of numerical solutions can not be compared with 
analytical ones though the former ones are very accurate. As an alternative to the numerical 
ones, semi-analytic expressions |^ and construction of certain RG invariant forms are useful for 
phenomenological analysis of the MSSM and its extensions 

The possibility of non-universality specific to Higgs masses was studied in a series of papers 
|[lO|,P], by noting constraints from b — > 57, cosmology and anomalous magnetic moment of muon 
and it was stressed that relaxing the scalar-mass universality assumption for the MSSM Higgs 
multiplets opens up many phenomenological possibilities (see also for Bs /i^ /i~ and cold 
dark matter issues related with the NUHM) . One of the aims of the present work is to present the full 
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form of semi-analytical expressions explicitly, so that all weak scale observables can be expressed in 
terms of GUT inputs. The analytical form of the results can provide considerable insight for similar 
issues (we ignore CP-violation during the numerical analysis of the NUHM, however, the full form 
of our results cover this issue too). Indeed, due to the complicated structure of the renormalization 
group equations , it is appealing to handle issues analytically and the solutions presented in this 
work can be useful for such an analysis even if they are given to the one loop order. As we will 
see, to keep the analysis simple, there are certain ignorance made on most of the correction terms, 
however, in the low tan (3 regime they do not affect our conclusions sizably and can further be added 
on demand. 

The outline of the rest of this work is as follows: In Section 2, we introduce our notation 
and conventions. In Section 3 we present the effects of non-universal Higgs masses terms on the 
supersymmetric mass spectra for varying tan (3 and scale values. A subsection of the same section 
is given to benchmark the semi- analytic results. Section 4 is devoted to our conclusions. The full 
form of the solutions of the RGEs can be found in the Appendix H. 



2. Notation and Conventions 

We define the basic parameters of the model as soft supersymmetry breaking scalar masses mp, 
gaugino masses M, the trilinear couplings Aq, bilinear coupling Bp and supersymmetric Higgs mass 
parameter /j,o, at the GUT scale. We assume third family dominance model and solve RGEs at the 
one- loop order. In this effective approach, by solving the RGEs explicitly, weak scale predictions are 
expressed in terms of GUT boundaries. We express Bino, Wino and Gluino with A/1.2,3, respectively, 
with a common initial value M. By writing the GUT boundaries, 

Ai = CAi Aq, Mj = CM J M, ruk ^ Ckmo (2.1) 

where i — t,b,T and j ~ 1, 2, 3, and for (mass)^ terms k = Hu, Hd, ii, tji, bn, tl, tr. We will express 
weak scale value of each quantity in terms of corresponding mSUGRA parameters toq, Aq, and 
a positive /i to be determined by the electroweak breaking conditions. From the solutions of the 
RGEs, weak scale and GUT scale values are connected and the most important restriction, in this 
respect, is the mass of Z boson: 

1 o o "Iff. — tan^ 3 rriif 

2^^^-^^+ \an^/3-l +^ (2.2) 

where tan/3 is the ratio of vacuum expectation values (vu/v^), and A stands for corrections on 
Higgs masses. We are interested in low vacuum expectation value (tan/3 = 10), for which complete 
list of semi- analytic solutions are given in the Appendix In addition to this, we will present 
graphical solutions for different tan /3 = 10 values in the next section. Instead of purely numerical 
values, expressing weak scale predictions in terms of GUT inputs proves very useful and helps to 
differentiate the importance of each term. In order to show the relative weigh of each term, we will 
express evolution of any soft (mass)^ term as in the following forms 

(mass)^ = 71 A^ -I- 72 Ao M + 73 + 74 c^^ m^ + 75 c|^ mg + 76 . (2.3) 

This decomposition enables one to lay stress upon the effects of non-universal Higgs mass choices. 
As it can be extracted from the above equation, sensitivity of each term to the initial values of mn^ 



and mn^ will be different. Notice that, by using the Mz constraint given in (2^) one can obtain n 
and this can be expressed as 



2|/^l' + m|,,, +m|,^ ^^^^ 



tan /3 -I- cot (3 
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hence, tree level relations of mass of physical Higgs boson can be written as in the foUowings 

= 2b/iiin2l3 (2.5) 
m^± = m\_^ + TO^ (2.6) 



m| ± 



\J (™Ao + "4)^ - 4m^^m| cos 2/? 



(2.7) 



Those relations will be modified, largely, due to top-stop loop corrections and /iq is the most affected 
one. Indeed, since the mass of the lightest CP-even Higgs boson is larger than 114 GeV ||l^ this 
correction must be included in the analysis. We will consider this correction and omit others in our 
effective approach. Meanwhile, the price that should be paid for that aim is predicting the spectra 
with small certain errors as will be shown in the following section. But this does not affect our 
conclusions, since the reaction of the SUSY particles to the non-universal Higgs boundary conditions 
is important for the present study. The necessary expression for the most important correction is 

^ ,,2 „,4 „• 4 



^«) -J^2^'yt In ( ^t:^ ) , (2.8) 



where ^ can be extracted from the following mass matrix 

fm^ + mf -f (i - |s^)M|cos2/3 to^ (At - cot /?) 



(2.9) 



y {At - ^ cot P) rnf + mf^ + 5 M| cos 2/3 



This 2x2 matrix can easily be diagonalized to obtain eigenvalues of stop quark masses in terms 
of GUT inputs, similarly the same should be done for m?, m? using the solutions presented in the 
appendices in order to get the full sparticle spectrum as usual {i.e. see [^). Indeed, having such 
analytic expressions is very useful to visualize the ingredients of sparticles to indirectly probe the 
allowed range of non-universality of Higgs bosons. As an example, for (tan/3 = 10) 

^ii 2 = -0-0523 Al + 0.192 AoM + 3.74 + (0.642 - 0.0176 c%_^ - 0.161 cjjj ml + 

- 0.245 AffT^^ (2.10) 

where the exact expression of 57 is a quite lengthy function of all terms appearing in the first line of 
( |2.10|) . Notice that, it can be obtained using the full forms of the solutions given in the Apppendix. 
Now, let us make a simplifying assumption /^o ^ ™o ~ M ^ Aq and fht ^ 2 toq on the part of 
( 2.10|) to approximately predict the composition of stop masses 

I ^ ~ -0.052 + 0.19 yloAf + 3.8 Af^ + (0.64 - 0.018 c^^ - 0.16 c^J + 

0.25 Af|=F 3.45 . (2.11) 



mi 

tl,2 



Using this analytical expression, for instance, one can conclude that weigh of up Higgs fields is 
larger than weigh of down Higgs fields but their relative weigh is negligible compared to other soft 



mass terms. To be specific, we will consider the specific reference point SPSla' 1 14 in the numerical 
analysis to benchmark the solutions provided. However, even under the above rough approximation 
we found mj^ = 472 GeV and m^^ — 506 GeV, to be compared with the exact results. For the mass 
spectra of SUSY particles, effects of up Higgs field can be dominant, however, as we will see in the 
next section this can not be generalized to other sectors. 

3. Numarical Analysis 

In this section tan/3 and scale evolutions of (mass)^ terms will be presented. For this aim, solutions 
of RGEs are performed such that high scale is set equal to 1.9 x 10^^ GeV and the supersymnietry 



- 3 - 



breaking scale is chosen as 1 TeV. With this choices, unification of the gauge coupUngs is satisfied 
at the GUT scale as gi = .92 = .gs = 0.718 ± 0.001. 

One can obtain the solutions of RGEs for any tan (3 ~ 10. To be specific, for tan /3 = 10, mass of 
the heavy SM fermions fix the Yukawa couplings at the same scale as Yt = 0.551, Yt — 0.0547, Yr = 
0.0685. As a brief summary of the semi-analytic solutions, this specific choice of tan/? yields the 
foUowings equations 



,2 



-0.102 Al + 0.375 Aq M - 1.93 - 0.709 ml + 0.0331 c^^ + 0.612 c^^ 



2 



mij^ = -0.0107 + 0.0309 AqM + 0.413 - 0.0241 + 0.955 c^j^ + 0.0333 c|^^_ 

TO?^ = -0.0367 Al + 0.134 Ao M + 4.33 + 0.757 ml + 0.00768 c^^ ml - 0.129 c^^ ml 

m?^ = -0.068^o + 0.25AoM + 3.15M2 + 0.527m^-0.0429c^^mg-0.194c^^m^ (3.1) 

TO?^ = -0.00534 yl^ + 0.0192 A^ M + 4.67 Af ^ + 0.988 + 0.0149 c]j^ ml - 0.0211 c]j^ ml 

mj^ = -0.00271 Al + 0.00216 Aq M + 0.493 + 0.994 - 0.0353 cf^^ + 0.0325 cj^^ ml 

= -0.00542 Al + 0.00432 ^0 M + 0.143 Af ^ + 0.989 ml + 0.0595 ml - 0.065 c^^ ml . 



Notice that the analytical expressions given in (3.1) arc constrained forms of the solutions presented 



in the Appendix (here we set and Ci — > 1, except for ch^ and ch^)- And they can be used 

at SPSla' point [|l^. We will benchmark our solutions using this point in the following subsection. 

Different scale and tan/3 effects can be extracted from the following figures (Figs. |l|-0). In Fig. 
1^, we show tan (3 and scale dependencies of the composition of . Normally, mass of up Higgs 
fields get contributions from any of the 28 terms given in (A. 2). When we assume CP is conserved 
~^ 0) and accept universality is in charge (except for Higgs fields), mass of the up Higgs field 
can be decomposed in a neat form as in (2.3) 

= li"-^ Al + 7f Ao M + 7f + 7f cf^^, ml + 7f cf^^, ml + ml . (3.2) 

As can be seen from both panels of the first figure, largest contribution to mass of up Higgs field 
comes from Gaugino sector (dashed-blue curves). Contribution of down Higgs field to up Higgs 
field is negligible, in other words, deviation of down Higgs from the universal choice can not yield 
a detectable effect on up Higgs field. In all Figs. 0-0, sohd red (green) curves corresponds to 



("^ffu) '-"^ related (mass) terms, which are mj^^, mj^^ 



,TOr„ ,m 



2 



and finally mf^, respectively. In order to show the effects of scale variations fix tan/3 = 10 (right 



panels) and for varying tan/3 values scale is fixed around the weak scale (left panels). The Figs. |1| 
denote that the gauge/gaugino sector contributions to scalar mass sector evolution increases scalar 
mass parameters as we go to the weak scale. It is visible in Figs. ^ and that a strong reaction can 
be detected in the slepton sector to non- universal Higgs mass terms and this is true for any tan/3 
value. Notice that this can be expected for Higgs bosons too (see Figs. |l]-||). We observe from Figs. 
1^-^ that, scalar top quarks are sensitive to NUHM only for very small tan /3 values (~ 2 — 3). 

During the numerical analysis we observed that following the physical Higgs bosons (except the 
CP-even light Higgs boson), sleptons are very sensitive to NUHM terms. Hence, we present Fig. ^ 
to show a bird-eye picture of the reaction of stau mass eigenvalues to NUHM parameters. As can 
be inferred from the very figure, reaction of sparticles to the mentioned non-universality drifts the 
mass predictions, to some extend. This effect ranges from a few GeV to ~ 30 — 40 GeV for different 
sparticles and it can be detectable since the correct spectrum is well known for the MSSM. See Tab. 
H for the reaction of the particles of the MSSM to NUHM. 

3.1 Benchmark of the solutions 

The most practical solution in order to test the trustability of our results is to use certain benchmark 
points. Though a large set of benchmark points and parameter lines in the MSSM parameter space 
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10 20 30 40 50 60 -0.2 -0.15 -0.1 -0.05 



tan (5 t 

Figure 1: Evolution of contributions for rn\j^ with tan/3 (left panel) and scale (right panel) in NUHM. 
Scale is shown with the dimensionless quantity t such that t = denotes the GUT scale and t ~ —0.2 
corresponds to the Z scale, tan /3 varies from 2 to 60. In both of the panels solid red, green and blue lines 
correspond to 74 " , 75 " and 7g " . Dashed red, green and blue lines correspond to " i 72 " ^^nd 
73^"' as given in (3.1) 




10 20 30 40 50 60 -0.2 -0.15 -0.1 -0.05 



tan (3 t 
Figure 2: The same as Fig. 1 but for rnlj^ 



10 20 



tan /? t 
Figure 3: The same as Fig. 1 but for rn^j^ 

is established, we will use one the the most studied points (sec for Snowmass Points and Slopes). 
Since we ignored most of the corrections except that of on the mass of the lightest CP-even Higgs 
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10 20 30 40 50 60 



tan /3 



Figure 4: The same as Fig. 1 but for 




Figure 5: The same as Fig. 1 but for m^^ 



20 30 



50 60 



tan /3 



-0.2 -0.15 -0.1 -0.05 

t 



Figure 6: The same as Fig. 1 but for m; 



boson, a strict comparison with the state of art programs hke ISAJET or SOFTSUSY 
shoul not be expected (see also and the given web page for onhne comparison). Nevertheless, 
resulting error should not be too high and there should be a visible correlation. We observed this 
is indeed the case for our semi-analytic solutions. To be definite, if tan/3 = 10, M — 250 GeV, 
mo — 70 GeV, Ao = —300 GeV, j^Ty 1 (which is the SPSla' reference point) then at the weak 
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Figure 8: 3D plot of variation stau mass eigenvalues, with non-universal coefficients at the SPSfa' point. 
The case ch^ = ch^ = 1 brings m^j = fff.4GeV (left) and mf^ = 199.7 GeV (right). Stop sensitivities are 



^«uMM _ 8 8 GeV and 5^,'^"" ~ 2.1 GeV. 



scale (at 1 TeV) we end up with Table ^ 

Comparison of these results with the reference point denotes that the errors in predicting 
m/j, TOf 2 s-nd rrij^o are negligible. For other mass terms errors are somewhat large especially in 
predicting mass of the lightest neutralino m^o^ ; here absolute error is ~ 8% which could be reduced if 
calculation are performed at two loops, corrections are noticed for all terms etc. However, apparent 
correlation is sufficient for our aim since we are basically interested in the reaction of those particles 
to the non-universal choices of the Higgs masses. Of course, this is true as far as corrections do 
not alter the weight of ch^ and cr^ on the SUSY particles, which we assumed to be true since 
the emerging mass difference of the worst prediction is less than ~ 8%. Nevertheless, a numerical 
simulation including all families and known corrections would be more decisive, which is beyond 
the scope of this work. 



4. Conclusions 

Using the semi-analytic solutions presented in this work it is observed that deviation from the 
universahty assumption of the Higgs fields does not induce serious problems as in the case of other 
soft (mass)^ terms (especially if ch^, ^ ch^ ^ 1). This can be inferred from Tab |l| in which 
coefficients of up and down Higgs fields are varied from to 2 mp. For this range, a striking 
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Particle 


SPSla' [GeV] 


This Work [GeV] 


% Difference 


^NUHM [Q^y] 


h° 


116.0 


110.3 


4.91 


0.2 




425.0 


425.7 


-0.17 


35.1 


A° 


424.9 


425.3 


-0.09 


35.1 




432.7 


432.8 


-0.02 


34.5 


ii 


366.5 


374.3 


-2.13 


5.1 




585.5 


578.9 


1.13 


5.1 


h 


506.3 


502.9 


0.67 


2.4 




545.7 


530.7 


2.75 


0.9 


fi 


107.9 


111.4 


-3.24 


8.8 




194.9 


199.7 


-2.46 


2.1 


X? 


97.7 


105.3 


-7.78 


0.2 


xl 


183.9 


194.3 


-5.65 


1.2 


xl 


400.5 


400.6 


-0.02 


15.9 




413.9 


417.3 


-0.82 


14.5 


xt 


183.7 


193.7 


-5.44 


1.3 


xt 


415.4 


417.9 


-0.60 


14.6 


Vr 


170.5 


176.6 


-3.58 


3.8 



Table 1: 

Numerical values for the mass of some of the supersymmetric particles and Higgs bosons in the 
reference point SPSla' and their comparison with our semi-analytical results. The third column 
is obtained by {SPSla' — our results) x 100/SPSla'. The fourth column denotes the sensitivity of 
each particle to the NUHM model parameters. The difference between maximal and minimal mass 
values is obtained by varying ch^ and ch^ in the [0, 2] interval and the emerging difference is called 
as sensitivity (JNUHM^ each term. 

difference can be observed on the mass of certain supersymmetric particles (like sleptons) while the 
others are insensitive to the mentioned phenomena (like lightest neutralino). 

Expected discovery of low energy SUSY at the Large Hadron Collider (LHC) and the In- 
ternational Linear Collider (ILC) will require reconstruction of the supersymmetric theory 
parameters from the experimental data. This is necessary not only for the minimal model but 
also for NUHM, especially if experimental data signalling deviations from the minimal supergravity 
model (mSUGRA) |2^ occurs. For this aim, precise measurements of mass of the light stau rrir^^, 
which is probably among the first sparticles to be discovered due to lepton nature and a light mass 
very sensitive to non-universality of the Higgs bosons, will be very suitable to shoot the NUHM 
parameter space. 
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A. Explicit Solutions for low tan/5 

In this part we present explicit form of our semi-analytic solutions which are obtained by solving 
the RGEs explicity, to the one loop order. The Gut scale is Mqut — 1-9 x 10"'^^ GeV and tan/3 — 10. 
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At the GUT scale we found the following results for gauge and Yukawa couplings 

gi = 0.7179, 52 = 0.7187, = 0.7195 
Yt = 0.5510, n = 0.0547, = 0.0685. 

(A.1) 

For the weak scale scale (~ 1 TeV), our results for soft (mass)^ terms read 

mlj^ = 0.000619^^4^ -7.8 x 10-^^4^ - 0.103 4^ + 0.00473 c^^^ 

+ 0.206 cli^ - 1.94 c|f3 + 0.0331 cjj^ ml + 0.612 cjj^ ml - 0.0319 4^ 
+ 0.0325 cl^ ml - 0.0325 c^^ mg - 0.387 ^ ml - 0.29 c^^ ml 

- 0.00572 CMi cos $12 - 0.0252 cmi cm^ M"^ cos $13 

- 0.0000612 Aq CAt, CM-, M cos + 2.13 x 10"'^ Aq ca^ cm, M cos $1^ 
+ 0.0122 Aq CAt CM, M cos $it - 0.168 cm-^ cm^ cos $23 

- 0.000535 Aq CAt CM2 M cos $26 + 1-12 X 10"*^ Aq CA^ CM2 M cos $2r 

+ 0.0726 Aq CAt CM2 M COS $2t - 0.00215 ^0 ca^ cm^ M cos $36 

+ 3.48 X 10"^ Aq CA^ CM3 M COS $3r + 0.293 Aq CAt ^Mg M COS $3t 

- 1.54 X lO'^AlcAtCA^ COS +0.000285 AqCA;, CAt C0S$t6 

- 3.29 X 10-^2 CA,CA, cos$t^ , (A.2) 

m]j^ = -0.00992 Al - 0.00272 AI + 0.000286 AI c\^ + 0.0361 c\ 

+ 0.449 c|f^ - 0.0613 cif3 + 0.955 ml + 0.0333 ml + 0.0224 cg^ 

- 0.0353 4^ ml + 0.0298 c^^ mo + 0.0232 Cj^ mo - 0.0642 c^^ mo 

- 0.000383 CMi CM2 cos $12 - 0.000749 cm, cms cos $13 

+ 0.000538 Aqca, CMi M cos$i6 + 0.000586 ca, ca/, M cos^i,- 

- 0.0000797 Aqca, cmi M cos$if - 0.0097 cos $23 

+ 0.0064 Aq Ca, M cos $2b + 0.0016 Aq Ca^ CM2 M cos <I>2r 

- 0.000762 Aq CAt CM2 M cos $2t + 0.0258 Aq ca, cm^, M cos $36 

- 0.0000721 Aqca, CM3 M cos$3r - 0.00307 Aqca, cm^ M cos $3* 
+ 0.0000554 CAb CA, cos + 0.00159 CAi, CAt cos^tb 

- 4.46 X 10"*^ CA, CAt cos^ti- , (A.3) 

mf^ = -0.0031 A^c^^ +5.35 x lO^^^A^c^^ -0.0342^^4^ - 0.00678 4fi 

+ 0.372 cIj^ + 4.04 c^^^, + 0.00768 c^^ m^ _ o.l29 4, "^0 ^ 0.014 c^^ 
+ 0.0108 c^^ ml - 0.0108 c^^ + 0.868 4 "^0 - 0.0965 4 "^o 

- 0.00197 CMi CAf2 cos $12 - 0.00865 cmi cm^ cos $13 

+ 0.00016 Aq ca, CMi ^ COS $16 - 1.29 x 10"^ Aq ca^ Cm, M cos $1^ 
+ 0.00403 Aq CAt cmi M cos $u - 0.0592 cms cms cos $23 

+ 0.00196 Aq CAi CM2 M COS $26 - 6.44 X 10"^ Aq Ca^ Cm2 M cos $2r 

+ 0.0239 Aq CAt CM2 M COS $2t + 0.00789 ca^ CM3 M cos $36 

- 0.000017 Aq ca^ CM3 M COS $3,- + 0.0965 Ao ca, cms M cos $34 
+ 0.0000106 ^2 CAb CA, cos $6^ + 0.000627 A^CA, CAt cos^tb 

- 1.17 X 10-M^CA,CAt cos$t^ , (A.4) 
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= 0.000412^^4^ -5.2 x IQ-'^ Alc\^ - 0.0686^0^1 + 0.0443 c^^^M^ 

- 0.168 cIj^ + 3.41 - 0.0429 cjj^ ml - 0.194 c|^^ + 0.0438 cg^ 

- 0.0434 c^^ + 0.0434 4„ ml - 0.193 c^^ TOq + 0-676 (?^^ ml 

- 0.00381 CMi Af ^ cos $12 - 0.0168 cmi CM3 cos $13 

- 0.0000408 ^oCAfc cmi M cos$i6 + 1.42 x 10"'^ ^0 ca, cmi M cos^i,- 
+ 0.0081 ^0 CAt cmi cos $u - 0.112 CM3 Af^ cos $23 

- 0.000357 ^0 CAy, cmi M cos $26 + 7.45 x 10"'^ Aq ca^ cm^ M cos $2t 
+ 0.0484 ^0 CAt CM2 M cos $2* - 0.00144 ca,, M cos $3^ 

+ 2.32 x 10"^ Ao ca^ CM3 M cos $3^ + 0.195 Aq ca^ cms M cos $34 

- 1.03 X 10^^ AlcA^CA^ cos $6^ + 0.00019 CAb CAt COS ^tb 

- 2.19 X lO-"^ AlcA^CAt cos^tT , (A.5) 



rn?^ = -0.00662 Al + 0.0000112 .4^ + 0.000191 Al c\ + 0.0162 c^^^ 

- 0.00483 c^^ + 4.66 c^^^ Af^ + 0.0149 c\_^ ml - 0.0211 ml + 0.972 c^^ ml 
+ 0.0217 c^^ - 0.0217 c^^ rr?,^ - 0.0279 c?^ + 0.0439 ml 

- 0.000119cMiCAf2 cos $12 - 0.000501 CMiCMjA^^ cos $13 

+ 0.000361 Aq ca^ cmiM cos $i6 - 2.71 x 10"*^ Aq ca^ cmi M cos $1.^ 

- 0.0000533 Ao CAt CMi M cos^u - 0.00647 cmj cm^ cos $23 
+ 0.00427 Ao CA„ CM2 M cos $26 - 0.0000136 ^0 ca, cms M cos $2r 

- 0.000509 Ao CAt CM2 Af cos $2t + 0.0172 ca^ cms cos $3^ 

- 0.0000363 Ao ca, cm^ M cos($3r - 0.00205 Ao ca^ cms M cos $3* 
+ 0.0000222 CAfc CA, cos $6,- +0.00106 ^^CAfc CAt cos^tb 

- 2.11 X 10"^ AgcA, CAt cos^tT , (A.6) 



m?^ =0.000011^24^-0.00274^24^-3.11 X 10"^ ^2 4^ +0.0365 c^^^ 

+ 0.457 Af^ + 0.000035 cli, - 0-0353 4^ ™o + 0-0325 4„ rnl + 0.0325 c^^ 
+ 0.965 cl^ ml + 0.0297 c^^ mf, + 0.0325 4 "^o - 0.065 c^^ ml 

- 0.000204 CMi CM2 Aif2 cos $12 + 3. x 10"^ cmi cms A'/^ cos $13 

- 3.46 x 10"*^ Ao CA^ CMiM cos + 0.00059 ^0 ca, cmi Af cos $1,- 
+ 2.42 x 10"'^ Ao CAt CMi M cos $u + 0.0000132 cm2 CM3 Af^ cos $23 

- 0.000014 ^0 CAfc CM2 M cos $26 + 0.00162 Ao ca^ cm^ M cos $2r 

+ 1.07 x 10"^ Ao CAt CM2 Af cos $2t - 0.0000176 ca^ cms Af cos $3^ 

- 0.0000178 ^oCA, CM3 Af cos$3t- + 1.78 x 10"® CAt cms Af cos $34 
+ 0.0000222 ^2 CA, CA, cos$6i--l-28 x 10"*^ ca, CAt cos$t6 

- 1.29 x IQ-'^AlcA^ CAt cos$tr , (A.7) 
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m?^ = 0.0000221 v4g 4^ ^ 0.00548 4^ -6.21 x lO-'^A^c^^ + 0. 147 c^^M^ 

- 0.00366 c^^ + 0.0000699 c^^ + 0.0595 cjj^ ml - 0.065 - 0.065 cl^ ml 
+ 0.0595 4^ ml + 0.929 c^^ ml - 0.065 c^^ + 0.13 c?^ TOq 

- 0.000409 cmi cos $12 + 5.99 x 10"*^ cmi cos $13 

- 6.93 X 10"^ Aq cmi M cos + 0.00118 ^0 ca, cmi M cos 

+ 4.83 X 10"'^ Ao cmi M cos $it + 0.0000264 cm^ cm^ M'^ cos $23 

- 0.000028 ^0 c^i, cm^M cos $26 + 0.00324 ^0 ca^ cm^M cos $2r 

+ 2.13 X 10"^ Ao CM2M cos $2t - 0.0000352 Aq ca^ cm^M cos $36 

- 0.0000355 Ao c^^ CM3 M cos $3^ + 3.57 x 10"^ Ao ca^ cms M cos $34 
+ 0.0000444 AgC^t Cyi, cos $6^-2.55 x lO"*^ Aq c^^ ca^ cos^tb 

- 2.58 X 10-'^ AlcA^CAt cos<^tT ■ (A.8) 
For Gauginos we found the foUowings 

Ml = 0.432 CMi Af, M2 = 0.833 Af, M3 = 2.51 cms M, . (A.9) 

Similarly, for trilinear terms we have 

At = -0.00198 AoCA, +3.81 x 10"^ Ao c^, + 0.27 Aq ca^ - 0.0303 cmi M - 0.231 cm^ M 
- 1.55 CM3M (A.IO) 

Ab = 0.147 AoCA, - 0.00041 Ao CA, - 0.0175 AoCa^ - 0.00484 cmi M - 0.0675 cms 

- 0.372 CM3 (A.ll) 

Ar = -0.00101 AoCAt +0.0989 Ac CA, + 0.0000811 AoCa^ - 0.0153 cmi M - 0.0493 M 
+ 0.00131 CM3M (A.12) 

Our expression for B is 

B = Bq- 0.0095 Ao CA, - 0.00276 Ao ca, - 0.354 Ao ca^ - 0.0301 cmi M - 0.371 M 
+ 0.518 CM3M (A.13) 

and, lastly, for the /j. parameter our result reads 

At = 0.995 MO- (A.14) 
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